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The mechanism of adsorption of the 1,4,5,8-naphthalene-tetracarboxylic-dianhydride (NTCDA) 
molecule on the Ag(110) surface is elucidated on the basis of extensive density functional theory 
calculations. This molecule, together with its perylene counterpart, PTCDA, are archetype organic 
semiconductors investigated experimentally over the past 20 years. We find that the bonding of the 
molecule to the substrate is highly site-selective, being determined by electron transfer to the LUMO 
of the molecule and local electrostatic attraction between negatively charged carboxyl oxygens and 
positively charged silver atoms in [110] atomic rows. The adsorption energy in the most stable site 
is 0.9eV. A similar mechanism is expected to govern the adsorption of PTCDA on Ag(110) as well. 

PACS numbers: 68.43.-h,73.20.Hb, 71.15Mb 



I. INTRODUCTION 

The adsorption of functional organic molecules on 
metal surfaces is of wide current interest both from tech- 
nological and fundamental points of view. Precise con- 
trol of the first monolayer is a prerequisite for fabricat- 
ing high-quality organic thin films, which are the basic 
constituents of hybrid organic-inorganic optical and elec- 
tronic devices^. Understanding molecule-molecule and 
molecule-substrate interactions is necessary to explain 
different surface architectures^, and knowledge of the 
local electronic structure^ and of mechanical properties 
is crucial for designing and controlling functionalities at 
the nanometer scale, for example in the fields of molecu- 
lar electronics 5,6 and molecular machines 7 . 




FIG. 1: (Color online) Local adsorption geometries studied: 
T (top), SB (short bridge), LB (long bridge) and H (hollow). 
The names refer to the position of the midpoint of the central 
C=C bond with respect to the underlying Ag(llO) surface. 



Among the many molecules whose interaction with 
solid substrates has been studied, NTCDA and its pery- 
lene counterpart, PTCDA, deserve special attention due 
to their unique properties. The aromatic cores of both 
molecules are terminated at opposite ends by two anhy- 
dride groups, as seen in Fig. ^ f° r NTCDA, and both 
of them have a D2h symmetry. Upon room tempera- 
ture deposition on silver surfaces they form commensu- 
rate superstructures, which depend on the orientation of 
the substrate^. Earlier NEXAFS studiea 10 i n i 12 of the 
adsorbed monolayers indicated that both molecules are 
parallel to the substrate and exhibit substrate-dependent 
spectral changes which were interpreted as evidence for 
covalent 7r-bonding\ The finding that such rather large 
molecules, which cover ~10-15 substrate atoms, can lock 
into preferred adsorption sites, had no simple explana- 
tion for a long time. Chemisorption via 7r-states delo- 
calized over the aromatic core, claimed to favor the high 
surface mobility required for self-assembly^, cannot eas- 
ily account for site-selectivity. Recently, Eremtchenko et 
a]±Lll proposed an interesting experimental argument 
that PTCDA binds to Ag(lll) via its central aromatic 
ring, but that this potential reaction center is inactive 
in the case of unsubstituted perylene. Later Hauschild 
et aim reported synchrotron-quality normal incidence 
X-ray standing wave (NIXSW) measurements and den- 
sity functional calculations on the same system, which 
revealed that anhydride side groups are slightly closer 
to the substrate than the aromatic core and also con- 
tribute to the bonding. Comparison of vibration spec- 
tra of adsorbed monolayers of PTCDA on different silver 
substrates by Tautz et a/ii 6 ^ 7 ^ showed, however, that 
beside common features there are also notable differences 
in local adsorption properties on different surfaces. The 
Ag(110) surface consists of a regular array of atomic rows 
and grooves in the [110] direction, and one expects (a) the 
lateral variation of the adsorption energy to be larger for 
the (110) surface than for (001) and (111) surfaces and 
(b) charge transfer to be more pronounced (see below). 
Similar differences should also exist for the adsorption of 
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NTCDA on silver substrates even though only the elec- 
tronic structure of the NTCDA/Ag(lll) interface has so 
far been studied in detailifi^. 

Our work is directly motivated by the experiments 
of Fink et al& who investigated the lateral ordering 
of NTCDA monolayers on several noble metal surfaces 
by means of LEED, STM and temperature-programmed 
desorption spectroscopy. Whereas two lateral superstruc- 
tures were found for NTCDA on Ag(lll) and Ag(OOl), 
only one simpler superstructure was observed for Ag(llO) 
and Cu(OOl). The NTCDA/Ag(110) interface, having a 
well-defined geometry and a relatively small number of 
atoms per unit cell, is therefore appealing for computa- 
tional work. 

Density functional theory (DFT) calculations 2 ^ have 
proven quite successful in describing covalent and ionic 
bonding of small molecules at surfaces^. However, it 
is not obvious whether bonding of large closed-shell or- 
ganic adsorbates on noble metals, where van-der-Waals 
interactions can play a significant role, will be adequately 
grasped by current semilocal implementations of DFT, 
like LDA and GGASS,. In general, since dispersion in- 
teractions are not taken into account when construct- 
ing the current density functionals^Si, one cannot ex- 
pect a good performance of DFT for pure van-der-Waals- 
bonded systems. It is unclear, however, how much van- 
der-Waals character do the molecule-metal bonds possess 
in the systems that concern us here. Even though the first 
DFT study of a large organic molecule interacting with 
a metal appeared almost ten years ago^£, still few such 
works exist. In the case of molecules like Cgo, which bond 
strongly to some noble metal substrates 26 , or when a sin- 
gle head group is involved in the bonding, as for alkane- 
thiols on Au(lll) 2 i, or when reactive metals such as N: 28 
are used as substrates, present-day DFT is largely suc- 
cessful. However, a recent study 2 ^ showed that calculated 
DFT adsorption energies of inert hydrocarbons are much 
smaller than experimental ones. For PTCDA, sometimes 
no bonding is found, even though experiments suggest a 
weak chemisorption 30,31 . On the more open and reactive 
Ag(llO) surface the molecule-metal bond may exhibit a 
more covalent (or ionic) character for which current den- 
sity functionals perform rather well. Fortunately, we do 
find that in our case DFT is successful, because much of 
the relevant experimental data, some extrapolated from 
related systems, is reproduced by calculations. 



II. COMPUTATIONAL DETAILS 

Calculations were performed with the ABINIT code 3 ^, 
using the PBE-GGA exchange-correlation functional 33 . 
Light atoms were described with Goedecker-Teter-Huttcr 
pseudopotentialsM, and silver atoms were treated within 
the Troullier-Martins scheme 3 ^ (s-component taken as lo- 
cal to avoid ghost states). The simulated system con- 
sisted of a NTCDA molecule on one side of a six-layer 
silver slab and a vacuum region equivalent to eight sil- 



ver layers. Thus the total computational box included 78 
atoms and a pool of 690 electrons, making it already a 
quite demanding system. Wave functions were expanded 
in plane waves with a 50Ry kinetic energy cutoff (60 and 
70Ry for test calculations - the results changed little) The 
Brillouin zone of the 3x3 surface unit cell was sampled by 
six special k-points 3 ^ together with a Fermi broadening 
of O.OlHa. 

The bottom three silver layers were fixed in their cal- 
culated bulk positions (PBE lattice constant 4.20 A), and 
the top three silver layers and the molecular degrees of 
freedom were relaxed according to Hellmann-Feynman 
forces, preserving the initial local symmetry. Geome- 
try optimization was performed using damped molecu- 
lar dynamics until the maximum force was below 10~ 3 
Ha/bohr. The potential energy surface describing the 
metal-molecule interaction is rather shallow, meaning 
that a substantial change in geometry leads only to a 
minor change in total energy. This motivated our strict 
convergence criterion. As test calculations, we calculated 
surface energies of clean silver surfaces, and obtained 
0.32, 0.41 and 0.62 eV per surface atom for (111), (001) 
and (110) faces in accord with previous GGA results and 
the experimental trend. The calculated work functions 
for these surfaces were 4.47, 4.27 and 4.13 eV, respec- 
tively, typical values for PBE using Troullier-Martins 
pseudopotentials. These values are also in accord with 
the experimental trend: the closed-packed (111) surface 
has the highest work function (4.46-4. 74eV) while the 
open (110) has the lowest work function among the low- 
index silver surfaces (4.14-4.52eV). This is quite impor- 
tant for the present system: the Ag(110) surface, having 
a smaller work function, can donate more charge to elec- 
tron acceptor molecules like NTCDA. These differences 
in work functions are certainly relevant for an interpreta- 
tion of adsorption properties of the molecules of interest 
on different substratesi 6 ^*^. 



III. RESULTS AND DISCUSSION 

After minimizing the residual forces, we obtained the 
following adsorption energies: 0.92eV for the T site, 
0.45eV for SB, and practically zero for LB and H. The 
adsorption energies were determined using the relation- 
ship E ad = -(E tot - E suh - E mo i), where E suh is the 
total energy of the bare slab (calculated using the same 
computational box), E mo \ is the total energy of an iso- 
lated NTCDA molecule (calculated using a large enough 
computational box), and E^ot is the total energy of 
substrate-adsorbate system. Positive adsorption energy 
means exothermic adsorption. The experimental adsorp- 
tion energy for NTCDA on Ag(110) is not known. For 
a superstructure of NTCDA on Ag(001) with a similar 

area per molecule (100. 2A 2 vs 106. 4A 2 ) a value of about 
l.OeV was obtained from the temperature-programmed 
desorption spectrum^, while for a denser superstructure 
of NTCDA on Ag(lll) an adsorption energy of l.leV was 




FIG. 2: Total electron density for the top site (upper row) and density differences for all four adsorption sites. Left column: 
zy-integrated densities (the units are eA ). Central column: averages over the [110] direction. Right: averages over the [001] 
direction. Solid lines represent electron charge accumulation regions, broken lines - depletion regions. 




FIG. 3: (Color online) Top: Contours of the total electron density (left) and of the electron density difference (right) in a plane 
0.5A below the NTCDA molecule at the top site; for An solid lines represent electron charge accumulation regions, broken 
lines - depletion regions. Bottom: 3D views of the LUMO (bi g symmetry), HOMO (i>3 9 ), and HOMO-1 (o u ) orbitals of the 
isolated NTCDA molecule. The distribution of accumulated charge closely resembles that of the LUMO of isolated NTCDA. 



deduced^. Thus, our value of 0.92eV is reasonable, but 
is most probably an underestimate, keeping in mind that 
usually GGA functional tend to underbind adsorbates 
on noble metals and that the adsorption energy on the 
Ag(llO) face should be even larger than those two val- 
ues. It is expected, however, that the differences between 
adsorption energies in different adsorption configurations 
should be more accurate. 

These differences are quite significant. The top site 
is predicted to be by far the most stable one. Fink et 
alA in fact also proposed that the top site, with all oxy- 
gens sitting almost on top of silver atoms underneath, 
should be the most favorable adsorption configuration 
(see Fig. 6 in Ref. 0). The proposal was motivated by 
analogy with the local adsorption geometry of PTCDA 
on Ag(llO), which had been determined in an elegant 
STM experiment, and noting that both PTCDA and 
NTCDA fit the Ag(llO) template (PTCDA spans four, 
and NTCDA spans three atomic rows - see Fig. ^ T site, 
and Fig. 4 in Ref. |37|). Our DFT calculations fully con- 
firm this hypothesis. For NTCDA on Ag(llO) the most 
favorable adsorption geometry is that for which oxygens 
lie directly above silver atoms in the [110] rows, just like 



for PTCDA on Ag(110). For the latter system an earlier 
computational study based on the empirical Universal 
Force Field (UFF)2£, predicted instead an equilibrium 
configuration equivalent to LB. The second most stable 
configuration, SB (Fig. 1), is also characterized by oxy- 
gen atoms on top of atomic rows but shifted half a lattice 
constant in the [110] direction. Although this remains to 
be explored in detail, configuration SB is likely a transi- 
tion state for diffusion along the Ag rows. In the case 
of the least favorable sites, LB and H, oxygen atoms 
sit between the close-packed rows. During the geome- 
try optimization, carboxyl oxygens moved towards the 
substrate, causing a bending of the molecule, just like 
for PTCDA on Ag(lll)i^. In the relaxed geometries, 
carboxyl oxygens were approximately 0.25 A and 0.30A 
closer to the silver substrate than the naphthalene core 
for the T and SB sites, respectively. For the T site, the 
carboxyl oxygens were situated 2.40A, the anhydride oxy- 
gens 2.50A, the carbons in the anhydride groups 2.54A, 
the carbons and hydrogens in the naphthalene core 2.62A 
from the topmost Ag(110) plane. An average distance of 
~ 3.0A was obtained for NTCDA on the less reactive 
Ag(lll) surface in a recent (NIXSW) experiment. On 
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Site 


T 


SB 


H 


LB 


Adsorption energy E a d [eV] 


0.92 


0.45 


~ 


~ 


Charge AN [electrons] 


0.38 


0.37 


0.34 


0.35 



TABLE I: Adsorption energies and net electron charges on 
the molecule determined from Eq. [5] after adsorption in four 
different sites. 

one hand, the distance between the NTCDA molecule 
and the more reactive Ag(110) surface should indeed be 
smaller than on the close-packed Ag(lll) surface. On 
the other hand, since the potential energy around the 
equilibrium molecule-metal distance is quite shallow in 
DFT, changes in that distance only lead to a tiny vari- 
ation of the adsorption energy. Thus we conclude that 
the calculated distances are subject to larger errors than 
the adsorption energies themselves. During the geom- 
etry optimization, in the case of the favorable sites (T 
and SB), silver atoms just below the carboxyl oxygens 
moved upwards by 0.07 A and 0.03 A with respect to the 
average level of the top layer, thus further reducing the 
corresponding Ag-0 distances. This naturally leads to 
the hypothesis that the interaction which leads to site- 
selectivity occurs between the carboxyl oxygens and the 
silver atoms in the [110] rows. If this is so, then what is 
the physical nature of this interaction? 

To answer this question, we analyze the electron den- 
sity redistribution caused by the adsorption of NTCDA 
at different sites by calculating the density difference 

An = n tot - n suh - n mo \ (1) 

where n tot is the total density of the system, and n sub and 
n mal are the densities of the substrate and molecule in 
their relaxed geometries. This function shows the charge 
rearrangement caused by adsorption of the molecule on 
the surface. Averages of the total densities and the den- 
sity differences along the [110] and [001] directions, as 
well as xy-integrated densities and density differences, 
are shown in Fig. [21 From these plots and 3D representa- 
tions of the density difference (not shown) one recognizes 
that (a) electron charge flows mainly from the metal to 
the molecule and (b) the charge accumulation region is 
closely related to the LUMO of the isolated molecule. 
This can clearly be seen in Fig. [21 where the total elec- 
tron density along with density difference in a plane 0.5A 
below the molecule (in the T configuration) is compared 
to the LUMO (big symmetry), HOMO (b3 g ) and HOMO- 
1 (a u ) orbitals of the free NTCDA molecule^. Prac- 
tically no back-donation from occupied orbitals occurs. 
For NTCDA on Ag(lll) NEXAFS studies of Gador et 
al. 1(111 showed a partial filling of the 7r-type LUMO. Our 
results are consistent with this finding. Various approx- 
imate charge partitioning schemes exist which associate 
charge with particular subsystems of the total system. 
We define the net electron charge on the molecule as 

AN = f An{z)dz, (2) 

J z 



where An(z) is the xy-integrated density difference, zq 
is the coordinate for which An(zo) is zero (see Fig. 
and z\ lies in the middle of the vacuum region, where all 
densities are practically zero. The net electron charge, as 
well as the adsorption energies for all four sites are sum- 
marized in Tabled Thus NTCDA acquires about 0.3-0.4 
electron charge upon adsorption. The charge transfer can 
be rationalized by the following simple argument. To a 
first approximation, the overall charge transfer is pro- 
portional to the difference of the chemical potentials of 
the substrate (referred to the vacuum level close to the 
surface of the metal) and of the molecule. The former 
is Msub = — "I'sub, ^sub being the work function of the 
Ag(110) surface, namely 4.13 eV (computed value). The 
latter, the so-called charge neutrality level^i, is approxi- 
mately equal to the midgap position^: 

&HOMO + clumo IP + EA 

Mmol ~ ^ ~ 2 ' 

where e#oMO> glumo are computed Kohn-Sham eigen- 
values while IP and EA are the ionization potential and 
the electron affinity of the isolated molecule. According 
to the last expression on the right hand side, fi mo i is just 
the negative of the absolute electronegativity defined by 
Mullikcn 43 . In our case, //ntcda ~ — 5.9eV, which agrees 
very well with the experimental value — 6.0eV—. Inter- 
estingly, the so-calculated absolute electronegativities are 
almost equal for NTCDA and PTCDA, even though the 
band gap is smaller for PTCDA than NTCDA. Thus, 
A*sub > Mntcda, so that electrons are transferred to the 
molecule, mainly to its LUMO. Such an interpretation 
is in agreement with chemical studies, since NTCDA is 
known to be an electron accepton4^. 

Charge transfer from the metal to the molecule leads to 
an increase of the work function of the combined system. 
From the total charge density distribution we obtain an 
increase A$ = +1.0eV. This corresponds to an induced 
negative electric dipole moment per NTCDA of 0.6eA 
which is slightly smaller than the crude estimate 0.4e x 
2.5 A partially because the relaxation of the ion cores 
makes an opposite contribution to the net dipole. No 
measurements of the work function for NTCDA/Ag(110) 
are available at present. Hill et al. 4 - found that the work 
function increases by +0.2eV upon adsorption of a re- 
lated molecule, 3,4,9, 10-perylenetetracarboxylic bisimi- 
dazole on polycrystalline Ag. A similarly small value was 
found for NTCDA on Ag(lll)±ML. Since the Ag(lll) 
surface has a higher work function, charge donation can 
be more substantial and A$ larger on the Ag(110) sur- 
face. However, our value of +1.0eV is probably over- 
estimated for two reasons. First, the calculated work 
function of the Ag(110) surface is slightly lower than the 
experimental one, and this favors charge donation to the 
LUMO of NTCDA. Second, as is well known, the band 
gaps of semiconductors^ and of organic molecules^ are 
underestimated in both LDA and GGA. The Kohn-Sham 
LUMO eigenvalue is usually lower in energy than the cor- 
responding electron affinity level and this also leads to a 
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a) Isolated subsytems b) Charge transfer to the LUMO 




d) Bending of the molecule c) Attraction between Ag and O 



FIG. 4: (Color online) Sketch of the mechanism of adsorp- 
tion of NTCDA on Ag(llO), viewed along the [110] direction, 
(a) The two isolated subsystems: note the electron excess on 
the peripheral oxygens and the positive charge on the naph- 
thalene core, (b) Charge transfer from the substrate to the 
molecule, mainly to its LUMO. (c) This transfer causes pos- 
itive charging of the topmost silver layer and a local electro- 
static attraction between silver atoms in the [110] rows below 
the negatively charged carboxyl oxygens, (d) The Ag-O at- 
traction causes a distortion of NTCDA from planarity, with 
peripheral oxygens moving closer to the Ag rows than the 
naphthalene core. 



stronger charge donation to the LUMO and thus to a 
larger A<£>. To cure these problems, expensive calcula- 
tions at the level of self-consistent many-body perturba- 
tion theory, e.g. the GW-approximation 5 — , must be per- 
formed. We believe, however, that the above-mentioned 
issues will not affect the proposed model for the site- 
selective adsorption of NTCDA on the Ag(110) surface. 

Table [I] shows that the overall charge on the molecule 
is very similar for all the adsorption sites studied and 
does not reflect the significant differences in adsorption 
energies. However, one notices from Fig. |2 that in the 
case of the favorable T and SB sites the charge trans- 
fer is stronger and more localized between the carboxyl 
oxygens and silver atoms underneath. Furthermore, the 
electrostatic energy gain crudely estimated by treating 
the molecule-metal system as a plane-parallel capacitor 
amounts to only 0.12 eV, i.e. much less than the ac- 
tual electrostatic contribution to the adsorption energy 
in those favorable sites. This also speaks for the fact that 
site selectivity is determined by the proximity of carboxyl 
oxygens and silver atoms. 

Our results suggest the following adsorption scenario 
of NTCDA on Ag(110) surface which leads to site-specific 



bonding. As a result of charge transfer to the LUMO of 
the molecule, silver atoms in the topmost layer become 
slightly positively charged. Carboxyl oxygens, being neg- 
atively charged even in the free molecule (~ 0.2 electron 
excess), are then attracted by electrostatic forces to the 
substrate. This attraction is maximal when the Ag-0 dis- 
tance is minimal, that is in the T and SB configurations. 
Thus, changes in the electron density are bigger and more 
localized for these two sites (Fig. 01 , as compared to the 
unfavorable LB and H configurations. The local Ag-0 at- 
traction leads to a distortion of NTCDA from planarity, 
i.e., the carboxyl oxygens move closer to the substrate 
than the naphthalene core. Most of these effects, to- 
gether with changes in bond lengths which reflect the 
density distribution of the LUMO, have been found and 
discussed for PTCDA on Ag(lll)A£. It is interesting to 
note the qualitative similarity between those two systems 
despite the fact that both molecules and substrates are 
different. Bonds through which the nodal planes of the 
LUMO pass are slightly elongated (maximum 0.03Afor 
the C=0 bond), and the others are slightly contracted. 
However, as mentioned in the introduction, the bond- 
ing is stronger and more site specific for the more open 
Ag(110) surface. Furthermore, free electron-like surface 
states of the Ag(lll) surface might also play a role in 
adsorption on that substrate. Figure 0] summarizes the 
mechanism of the bonding and bending of the molecule 
on the Ag(110) surface. 

Concerning the lateral ordering of the molecules, we 
studied adsorption in the favored T configuration for two 
different lateral surface unit cells: (3,0/0,3) and the ex- 
perimentally found (3,0/1,3)2.. Both unit cells contain 
one molecule and the same number of substrate atoms 
and both structures are rather open. We found that 
the adsorption energies in both of these arrangements 
are essentially the same. This happens, most probably, 
because semilocal GGA functionals fail to describe long- 
range van-der-Waals interactions^, which are likely to be 
most important for the preference of one superstructure 
over the other. 



IV. CONCLUSIONS 

In conclusion, we have studied site-selective adsorption 
properties of NTCDA on the Ag(110) surface by means 
of density functional calculations. The interaction of the 
molecule with the substrate can be summarized by the 
following scenario: (1) NTCDA, being an electron accep- 
tor, takes ~ 0.4 electron from the silver substrate and 
the LUMO is partially filled; (2) this transfer leads to a 
slight positive charging of the silver atoms in the topmost 
layer and weakening of the C=0 bonds, thus enhanc- 
ing the local electrostatic interaction between Ag atoms 
in the Ag[110] rows and carboxyl oxygens (this is the 
main difference between favorable and unfavorable sites); 
(3) the Ag-0 attraction distorts the molecule from the 
planar configuration. Being quite similar to the mecha- 
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nism recently established for PTCDA on Ag(lll)i^, we 
believe that the proposed scenario is also applicable to 
the adsorption of PTCDA on Ag(110)2£. Experimentally 
found differences between different substratesiSiiLiS are 
still to be explored in detail. Very recently, a single po- 
lar N-Cu bond was found to be responsible for bonding 
of inclined adenine on Cu(110)5i, like for alkane-thiols 
on Au(lll)2i, suggesting that local electrostatic bonding 
can occur for many organic molecules on open fcc(llO) 
surfaces. In the case of aromatic acceptor molecules with 
several peripheral electronegative atoms (O, N, halogens) 
adsorbed flat on noble metal fcc(llO) surfaces, the origin 
of site-selectivity is less obvious. Our results suggest that 
observed commensurate superstructures can arise from 
the close match of those atoms and electron-depleted 



substrate atoms. DFT computations appear capable of 
grasping lateral variations of the resulting electrostatic 
interactions. 
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